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Abstract. The formation of the first galaxies at redshifts z ~ 10 — 15 signaled the transition from the simple initial state of 
the universe to one of ever increasing complexity. We here review recent progress in understanding their assembly process 
with numerical simulations, starting with cosmological initial conditions and modelling the detailed physics of star formation. 
In particular, we study the role of HD cooling in ionized primordial gas, the impact of UV radiation produced by the first 
stars, and the propagation of the supernova blast waves triggered at the end of their brief lives. We conclude by discussing 
promising observational diagnostics that will allow us to probe the properties of the first galaxies, such as their contribution 
to reionization and the chemical abundance pattern observed in extremely low-metallicity stars. 
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INTRODUCTION 

One of the key goals in modern cosmology is to study 
the assembly process of the first galaxies, and under- 
stand how the first stars, stellar systems, and massive 
black holes formed at the end of the cosmic dark ages, 
a few hundred million years after the Big Bang. With the 
formation of the first stars, the so-called Population III 
(Pop III), the universe was rapidly transformed into an 
increasingly complex, hierarchical system, due to the en- 
ergy and heavy element input from the first stars and ac- 
creting black holes [7, 17, 26, 64]. Currently, we can di- 
rectly probe the state of the universe roughly a million 
years after the Big Bang by detecting the anisotropics in 
the cosmic microwave background (CMB), thus provid- 
ing us with the initial conditions for subsequent struc- 
ture formation. Complementary to the CMB observa- 
tions, we can probe cosmic history all the way from the 
present-day universe to roughly a billion years after the 
Big Bang, using the best available ground- and space- 
based telescopes. In between lies the remaining frontier, 
and the first galaxies are the sign-posts of this early, for- 
mative epoch. 

The tight correlation between the central black hole 
mass and properties of the host galactic spheroid ob- 
served at low redshift indicates that galaxy formation 
and massive black hole growth are fundamentally re- 
lated [35, 36, 39, 96]. A popular class of models to ex- 
plain this coiTelation invokes black hole feedback in the 
form of winds or thermal energy input to regulate gas 
fueling [86]. Such feedback effects might be even more 
important in the shallow potential wells of high-redshift 
galaxies [57] . An open problem is to understand the early 
growth rate of black holes, before they reach 10^ M©, 



as a function of their cosmic environment, a process that 
ultimately leads to the formation of quasar-sized massive 
black holes at z > 6 [57] . Due to the complexity involved, 
one needs to carry out high-resolution numerical simula- 
tions to investigate both the build-up of the first stellar 
systems and the early growth of massive black holes. 

The internal structure of the first galaxies involves a 
wide range of physical scales, from the virial radius of 
the host dark matter halo to the much smaller scales such 
as the Jeans length of the cooling gas and the Bondi 
radius of the first massive black hole. Furthermore, the 
dynamical evolution of the first galaxies will involve 
turbulence and shocks that must be resolved. Some of 
the key questions are: What is the angular momentum 
content of the cold gas, and what processes affect angular 
momentum transport in the first rotationally-supported 
gas flow? Do global gravitational instabilities, such as 
bars-within-bars, give rise to rapid transport, and what 
are the consequences of rapid transport for the growth of 
the central object? What is the role of local gravitational 
instability and fragmentation? What is the impact of 
the radiation produced by the accreting black hole on 
the thermodynamic and chemical evolution of the first 
galaxy? At which point will a multi-phase interstellar 
medium form? 

To simulate the build-up of the first stellar systems, 
we have to address the feedback from the very first stars 
on the surrounding intergalactic medium (IGM), and the 
formation of the second generation of stars out of ma- 
terial that was influenced by this feedback. There are a 
number of reasons why addressing the feedback from the 
first stars and understanding second-generation star for- 
mation is crucial: 

(i) The first steps in the hierarchical build-up of struc- 



ture provide us with a simplified laboratory for studying 
galaxy formation, which is one of the main outstanding 
problems in cosmology, given that we have successfully 
determined the parameters of the expanding background 
universe. The overall strategy is to start with cosmologi- 
cal initial conditions, follow the evolution up to the for- 
mation of a small number (A'^ < 10) of Pop III stars, and 
trace the ensuing expansion of the supernova (SN) blast 
waves after they died together with the dispersal and 
mixing of the first heavy elements, towards the forma- 
tion of second-generation stars out of enriched material 
[42, 101]. 

( ii) The initial burst of Pop III star formation may have 
been rather brief due to the strong negative feedback ef- 
fects that likely acted to self-limit this formation mode 
[41, 104]. Second-generation star formation, therefore, 
might well have been cosmologically dominant com- 
pared to Pop III stars. Despite their importance for cos- 
mic evolution, e.g., by possibly constituting the majority 
of sources for the initial stages of reionization at z > 10, 
we currently do not know the properties, and most im- 
portantly the typical mass scale, of the second-generation 
stars that formed in the wake of the very first stars. 
(Hi) A subset of second-generation stars, those with 
masses below ~ 1 Mq, would have survived to the 
present day. Surveys of extremely metal-poor Galactic 
halo stars therefore provide an indirect window into the 
Pop III era by scrutinizing their chemical abundance pat- 
terns, which reflect the enrichment from a single, or at 
most a small multiple of. Pop III SNe [9, 37]. Stellar ar- 
chaeology thus provides unique empirical constraints for 
numerical simulations, from which one can derive theo- 
retical abundance patterns to be compared with the data. 

Recently, observations of the highest-redshift quasars 
and galaxies to date have opened up a tantalizing window 
into the state of the universe at z > 7 [8, 13, 32, 48, 55, 
76, 80, 92, 103]. Studying the nature of star formation 
at even higher redshifts is of great interest, because the 
systems that we can directly observe at z > 7 will already 
show the signature of previous stars. An intriguing exam- 
ple is the recently discovered 7— band dropout HUDF- 
JD2, found through deep HST/YW/Spitzer imaging, 
and interpreted to be a massive (~ 6 x lO'^ Mq) post- 
starburst galaxy at z > 6.5 [65]. Since the observed stel- 
lar population in HUDF-JD2 is old, star formation could 
have begun already at z ~ 15, and this galaxy could 
thus contain the signature of first- and second-generation 
stars. 

Existing and planned observatories, such as HST, 
Keck, VLT, and the James Webb Space Telescope 
(JWST), planned for launch around 2013, yield data 
on stars and quasars less than a biUion years after the 
Big Bang. The ongoing Swift gamma-ray burst (GRB) 
mission provides us with a possible window into massive 
star formation at the highest redshifts [18, 22, 56]. Mea- 



surements of the near-IR cosmic background radiation, 
both in terms of the spectral energy distribution and 
the angular fluctuations provide additional constraints 
on the overall energy production due to the first stars 
[29, 33, 52, 62, 79]. Understanding the formation of the 
first galaxies is thus of great interest to observational 
studies conducted both at high redshifts and in our local 
Galactic neighborhood. 

ASSEJVIBLY OF THE FIRST GALAXIES 

How massive were the first galaxies, and when did they 
emerge? Theory predicts that dark matter (DM) halos 
containing a mass of ~ 10^ Mq and collapsing at z ~ 
10 — 15 were the hosts for the first bona fide galaxies. 
These dwarf systems are special in that their associated 
virial temperature exceeds the threshold, ~ 10^ K, for 
cooling due to atomic hydrogen [70]. These so-called 
'atomic-cooling halos' did not rely on the presence of 
molecular hydrogen to enable cooling of the primordial 
gas. In addition, their potential wells were sufficiently 
deep to retain photoionization heated gas, as well as SN 
shocked gas, in contrast to the shallow potential wells of 
minihalos [28, 42, 61, 66]. These are arguably minimum 
requirements to set up a self -regulated process of star for- 
mation that comprises more than one generation of stars, 
and is embedded in a multi-phase interstellar medium. In 
the following, we discuss some of the key processes that 
govern the assembly of these first dwarf galaxies. 

To avoid confusion, we would hke to cormnent on 
a change in terminology. It has become evident that 
Pop III star formation might actually consist of two dis- 
tinct modes: one where the primordial gas collapses into 
a DM minihalo (see below), and one where the metal- 
free gas becomes significantly ionized prior to the onset 
of gravitational runaway collapse [49]. We had termed 
this latter mode of primordial star formation 'Pop II. 5' 
[41, 49, 60]. To more clearly indicate that both modes 
pertain to metal-free star formation, we here follow the 
new classification scheme suggested by Chris McKee 
(see McKee in these proceedings). Within this scheme, 
the minihalo Pop III mode is now termed Pop III.l, 
whereas the second mode (formerly 'Pop II. 5') is now 
called Pop III.2. The hope is that McKee's terminology 
wiU gain wide acceptance. 



Role of HD Cooling 

While the very first Pop III stars (so-called Pop III.l), 
with masses of the order of 100 Mq, formed within DM 
minihalos in which primordial gas cools by H2 molecules 
alone, the HD molecule can play an important role in 



the cooling of primordial gas in several situations, allow- 
ing the temperature to drop well below 200 K [1, 14]. In 
turn, this efficient cooling may lead to the formation of 
primordial stars with masses of the order of 10 M© (so- 
called Pop III. 2) [49]. In general, the formation of HD, 
and the concomitant cooling that it provides, is found 
to occur efficiently in primordial gas which is strongly 
ionized, owing ultimately to the high abundance of elec- 
trons which serve as catalyst for molecule formation in 
the early universe [84]. 

Efficient cooling by HD can be triggered within the 
relic H II regions that surround Pop III.l stars at the 
end of their brief lifetimes, owing to the high electron 
fraction that persists in the gas as it cools and recombines 
[5 1 , 67, 1 05] . The efficient formation of HD can also take 
place when the primordial gas is collisionally ionized, 
such as behind the shocks driven by the first SNe or in 
the viriaUzation of massive DM halos [41, 49, 59, 85]. 
In Figure 1, we show the HD fraction in primordial gas 
in four distinct situations: within a minihalo in which 
the gas is never strongly ionized, behind a 100 km s~^ 
shock wave driven by a SN, in the virialization of a 
3a DM halo at redshift z = 15, and in the reUc H II 
region generated by a Pop III. 1 star at z ~ 20 [49]. Also 
shown is the critical HD fraction necessary to allow the 
primordial gas to cool to the temperature floor set by the 
CMB at these redshifts. Except for the situation of the 
gas in the virtually un-ionized minihalo, the fraction of 
HD becomes large quickly enough to play an important 
role in the cooling of the gas, allowing the formation of 
Pop III.2 stars. 

Figure 2 schematically shows the characteristic 
masses of the various stellar populations that form in 
the early universe. In the wake of Pop III.l stars formed 
in DM minihalos. Pop III.2 star formation ensues in 
regions which have been previously ionized, typically 
associated with relic H II regions left over from massive 
Pop III. 1 stars collapsing to black holes, while even later, 
when the primordial gas is locally enriched with metals, 
Pop II stars begin to form [20, 41]. Recent simulations 
confirm this picture, as Pop III.2 star formation ensues 
in relic H II regions in weU under a Hubble time, while 
the formation of Pop 11 stars after the first SN explosions 
is delayed by more than a Hubble time [42, 105, 106]. 
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FIGURE 1. Evolution of the HD abundance, Xhd> in pri- 
mordial gas which cools in four distinct situations. The solid 
line corresponds to gas with an initial density of 100 cm^^, 
which is compressed and heated by a SN shock with velocity 
Vsh = 100 km s~' at z = 20. The dotted line corresponds to 
gas at an initial density of 0.1 cm~^ shocked during the for- 
mation of a 3<J halo at z = 15. The dashed line corresponds 
to unshocked, un-ionized primordial gas with an initial density 
of 0.3 cm"-' collapsing inside a minihalo at z = 20. Finally, 
the dash-dotted line shows the HD fraction in primordial gas 
collapsing from an initial density of 0.3 cm"^ inside a relic 
H II region at z = 20. The horizontal line at the top denotes 
the cosmic abundance of deuterium. Primordial gas with an 
HD abundance above the critical value, XuDjCrit' denoted by the 
bold dashed line, can cool to the CMB temperature witliin a 
Hubble time. 



■M 



3i» 



\ 

Pop. 111.2" 



S 



Radiative Feedback 

Due to their extreme mass scale. Pop III.l stars emit 
copious amounts of ionizing radiation, as well as a strong 
fiux of H2-dissociating Lyman- Werner (LW) radiation 
[16, 81]. Thus, the radiation from the first stars dramati- 
cally influences their surroundings, heating and ionizing 
the gas within a few kiloparsec (physical) around the pro- 



FIGURE 2. Characteristic stellar mass as a function of 
redshift. Pop III.l stars, formed from unshocked, un-ionized 
primordial gas are characterized by masses of the order of 
100 Mq . Pop II stars, formed in gas which is enriched with met- 
als, emerged at lower redshifts and have characteristic masses 
of the order of 1 M©. Pop III.2 stars, formed from ionized 
primordial gas, have characteristic masses reflecting the fact 
that they form from gas that has cooled to the temperature of 
the CMB. Thus, the characteristic mass of Pop III.2 stars is a 
function of redshift, but is typically of the order of 10 Mq. 



genitor, and destroying the H2 and HD molecules locally 
within somewhat larger regions [2, 5, 34, 51, 54, 98]. 
Additionally, the LW radiation emitted by the first stars 
could propagate across cosmological distances, allow- 
ing the build-up of a pervasive LW background radiation 
field [43]. 

The impact of radiation from the first stars on their lo- 
cal surroundings has important implications for the num- 
bers and types of Pop III stars that form. The photo- 
heating of gas in the minihalos hosting Pop III.l stars 
drives strong outflows, lowering the density of the pri- 
mordial gas and delaying subsequent star formation by 
up to 100 Myr [51, 98, 105]. Furthermore, neighboring 
minihalos may be photoevaporated, delaying star forma- 
tion in such systems as well [4, 42, 83, 93, 99]. The pho- 
todissociation of molecules by LW photons emitted from 
local star-forming regions will, in general, act to delay 
star formation by destroying the main coolants that al- 
low the gas to collapse and form stars. 

The photoionization of primordial gas, however, can 
ultimately lead to the production of copious amounts of 
molecules within the relic H II regions surrounding the 
remnants of Pop III.l stars [50, 67, 70, 75]. Recent simu- 
lations tracking the formation of, and radiative feedback 
from, individual Pop III. 1 stars in the early stages of the 
assembly of the first galaxies have demonstrated that the 
accumulation of relic H II regions has two important ef- 
fects. First, the HD abundance that develops in relic H II 
regions allows the primordial gas to re-collapse and cool 
to the temperature of the CMB, possibly leading to the 
formation of Pop III. 2 stars in these regions [51, 106]. 
Second, the molecule abundance in relic H II regions, 
along with their increasing volume-filling fraction, leads 
to a large optical depth to LW photons over physical dis- 
tances of the order of several kiloparsecs. The develop- 
ment of a high optical depth to LW photons over such 
short length-scales suggests that the optical depth to LW 
photons over cosmological scales may be very high, act- 
ing to suppress the build-up of a background LW radia- 
tion field, and mitigating negative feedback on star for- 
mation. 

Figure 3 shows the chemical composition of primor- 
dial gas in relic H II regions, in which the formation of 
H2 molecules is catalyzed by the high residual electron 
fraction. Figure 4 shows the average optical depth to LW 
photons across the simulation box, which rises with time 
owing to the increasing number of relic H II regions. 



IMassive Black Hole Growth 

The observations of quasars at redshift z > 6 in the 
Sloan Digital Sky Survey (SDSS) suggest that some black 
holes in the early universe grew to have masses in excess 




FIGURE 3. The chemical interplay in relic H II regions. 
While all molecules are destroyed in and around active H II 
regions, the high residual electron fraction in relic H II re- 
gions catalyzes the formation of an abundance of H2 and HD 
molecules. The light and dark shades of blue denote regions 
with a free electron fraction of 5 x lO^"' and 5 x 10^**, respec- 
tively, while the shades of green denote regions with an H2 
fraction of lO""*, 10"^ and 3 x 10"^, in order of decreasing 
brightness. The regions with the highest molecule abundances 
lie within relic H II regions, which thus play an important role 
in subsequent star formation, allowing molecules to become 
shielded from photodissociating radiation and altering the cool- 
ing properties of the primordial gas. 

of 10^ Mq within the first billion years after the Big 
Bang [30, 31]. Many massive Pop III stars may have 
collapsed directly to form black holes at the end of 
their brief lifetimes, thus providing the seed black holes 
which then accreted gas and grew to be the supermassive 
black holes powering the high-redshift quasars that are 
observed today [45]. 

In order for Pop III remnant black holes to accrete 
'^10^ M0 of mass by z ^ 6, these objects would have 
to accrete gas at or near the Eddington rate for hundreds 
of millions of years. In turn, these black holes can only 
accrete at such high rates if the surrounding density ex- 
ceeds 10^ cm^^ (see Figure 5) [50]. However, as Figure 6 
shows, the photoheating of the gas surrounding the first 
Pop III stars in minihalos drives strong outflows from the 
central star, leaving a remnant black hole to reside in the 
middle of the evacuated minihalo with densities typically 
well below 1 cm^^. Thus, the high accretion rates neces- 
sary to explain the presence of 10^ M0 black holes at 
z ~ 6 are not possible for at least 50 Myr after the forma- 
tion of the black hole, at which time DM halo mergers 
and the re-collapse of the photoheated gas may raise the 
central density of the gas to above 10^ cm^^ [51]. Sim- 
ulations tracking the later accretion history of such seed 
black holes have also found that it is difficult for these 
objects to grow fast enough to explain the SDSS quasar 
observations (see also Alvarez et al. in these proceed- 
ings) [57, 74]. Therefore, more exotic origins for these 
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FIGURE 4. Optical depth to LW photons averaged over two 
different scales, as a function of redshift. The diamonds denote 
the optical depth averaged over the entire cosmological box 
of comoving length 660 kpc, while the plus signs denote the 
optical depth averaged over a cube of 220 kpc (comoving) per 
side, centered on the middle of the box. The solid line denotes 
the average optical depth that would be expected for a constant 
H2 fraction of 2 x 10^^ (primordial gas), which changes only 
due to cosmic expansion. 
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FIGURE 6. The density, temperature and radial velocity of 
primordial gas heated and ionized by radiation from a 100 Mq 
and 200 Mq, primordial star after their main-sequence lifetimes 
of 3 Myr and 2 Myr, respectively. The central density drops to 
well below 1 cm^ , which is insuffiecient to enable subsequent 
accretion at the Eddington rate (see Figure 5). 
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FIGURE 5. The mass of an initially 100 M@ black hole as 
a function of time, assuming the black hole accretes gas at 
a temperature of 200 K and at constant density. The solid, 
dotted, dashed, and dot-dashed lines show cases with densities 
of 0.1 cm^-^, 1 cm^-^, 10 cm^^, and 100 cm^-', respectively. 
The triple dot-dashed line shows the mass of the black hole as 
a function of time, assuming that it accretes at the Eddington 
limit. Clearly, for the black hole to begin accreting at such a 
high rate while its mass is still of the order of 100 Mq, the 
surrounding density must exceed 10^ cm^^. 



supermassive black holes may be required [10, 19, 58]. 



The First Supernova Explosions 

Recent numerical simulations have indicated that pri- 
mordial stars forming in DM minihalos typically attain 
100 Mq by efficient accretion, and might even become as 
massive as 500 Mq [21,71, 72, 107]. After their main- 
sequence lifetimes of typically 2 — 3 Myr, stars with 
masses below 2± 100 Mq, are thought to collapse directly 
to black holes without significant metal ejection, while 
in the range 140 — 260 Mq a pair-instability supernova 
(PISN) disrupts the entire progenitor, with explosion en- 
ergies ranging from 10^' — 10^^ ergs, and yields up to 
0.5 [45, 46]. Less massive primordial stars with a high 
degree of angular momentum might explode with simi- 
lar energies as hypernovae [94, 97]. 

The significant mechanical and chemical feedback ef- 
fects exerted by such explosions have been investigated 
with a number of detailed calculations, but these were ei- 
ther performed in one dimension [53, 59, 77], or did not 
start from realistic initial conditions [23, 69]. The most 
realistic simulation to date took cosmological initial con- 
ditions into account, and followed the evolution of the 
gas until the formation of the first minihalo at z ~ 20 
(see Figure 7) [42]. After the gas approached the 'loiter- 
ing regime' at nn — lO'* cm^^, the formation of a pri- 
mordial star was assumed, and a photoheating and ray- 
tracing routine determined the size and structure of the 
resulting H II region (see Figure 8) [51]. An explosion 
energy of 10^^ ergs was then injected as thermal energy 
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FIGURE 7. Hydrogen number density averaged along the 
line of sight in a slice of 14 kpc (comoving) around the first star, 
forming in a halo of total mass Myj^ ~ 5 x 10^ Mtj at z ~ 20. 
Evidently, the host halo is part of a larger group of less massive 
minihalos, and subject to the typical bottom-up evolution of 
structure formation. 




Redshift: 

z = 21 

Boxsize: 
150/h kpc (comoving) 

Slice Width: 
10/h kpc (comoving) 



FIGURE 8. Temperature averaged along the line of sight in 
a slice of 14 kpc (comoving) around the star after its main 
sequence lifetime of 2 Myr. Ionizing radiation has penetrated 
nearby minihalos and extends up to 5 kpc (physical) around 
the source, heating the IGM to roughly 2 x 10^ K, while some 
high-density regions have effectively shielded themselves. 



into a small region around the progenitor, and the sub- 
sequent expansion of the SN remnant was followed until 
the blast wave effectively dissolved into the IGM. The 
cooling mechanisms responsible for radiating away the 
energy of the SN remnant, the temporal behavior of the 
shock, and its morphology could thus be investigated in 
great detail (see Figure 9). 



Dynamical Evolution 

In analogy to present-day SNe, the evolution of 
SN remnants in the early universe can be decom- 
posed into four physically distinct stages: free expan- 
sion (FE), Sedov-Taylor blast wave (ST), pressure-driven 
snowplow (PDS), and momentum-conserving snowplow 
(MCS) [73]. This allows the introduction of a simple an- 
alytic model, which reproduces the simulation results of 
[42] relatively well (see Figure 9). It was found that the 
SN remnant propagates for a Hubble time at z ~ 20 to a 
final mass-weighted mean shock radius of 2.5 kpc (phys- 
ical), roughly half the size of the H II region, and sweeps 
up a total (gas) mass of 2.5 x 10^ Mp,. The radial disper- 
sion of the SN remnant increased dramatically once the 
shock left the host halo, although the bulk of the swept- 
up gas was expelled into the voids of the general IGM. 



Mechanical Feedback 

As confirmed in [42], highly energetic SN explosions 
are sufficient to entirely disrupt the host halo and evac- 
uate their gaseous content [60]. This is directly related 
to their shallow potential wells, and was also found in 
previous calculations [23, 53]. Subsequent star formation 
does not ensue for at least a Hubble time at z ~ 20, and 
in the specific scenario considered, the first Pop II stars 
are expected to form out of enriched material at z ~ 10 
[42]. 

Additional simulations in the absence of a SN explo- 
sion were performed to investigate the effect of photo- 
heating and the impact of the SN shock on neighbour- 
ing minihalos (see Figure 10) [42]. For the case dis- 
cussed in this work, the SN remnant exerted positive me- 
chanical feedback on neighbouring minihalos by shock- 
compressing their cores, while photoheating marginally 
delayed star formation. Although a viable theoretical 
possibility, secondary star formation in the dense shell 
via gravitational fragmentation was not observed, pri- 
marily due to the previous photoheating by the progen- 
itor and the rapid adiabatic expansion of the post-shock 
gas [59, 60, 77]. 
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FIGURE 9. The mass-weighted mean shock radius (Panel 
(a)), and swept-up (gas) mass (Panel (b)) of a 10^^ ergs SN 
explosion in the high-redshift universe as a function of time 
(black dots), compared to the analytic model (dashed line). In 
the specific scenario considered, a final mass-weighted mean 
shock radius of 2.5 kpc (physical) and a final swept-up mass of 
2.5 X was found. The shaded region shows the radial 

dispersion of the shock, indicating that it increases dramatically 
once the SN remnant leaves the host halo and encounters the 
first neighbouring minihalos. 



2.0 
1.8 
1.6 

J, 1.2 
1.0 
0.8 
0.5 



holo (c) 



hcio (b) 



halo (o) 



100 

t,K [My] 



FIGURE 10. The collapse times and distances from the SN 
progenitor for all star-forming minihalos affected by the SN 
shock. The shades of the symbols indicate their affiliation, 
i.e. black, dark grey and light grey symbols represent the no- 
feedback, photoheating-only and main simulation runs, respec- 
tively, while the shapes of the symbols denote the individual 
halos. For orientation, the dashed line shows the mass-weighted 
mean shock radius at late times according to Figure 9. In the 
cosmological realization considered, photoheating significantly 
delays star formation, while the SN shock compresses gas in 
neighbouring minihalos and slightly accelerates their collapse. 



FIGURE 11. Hydrogen number density averaged along the 
line of sight in a slice of 14 kpc (comoving), overlayed with the 
distribution of all metal particles (bright orange) after 200 Myr, 
when the shock stalls. Most metals are dispersed into the voids 
of the general IGM, since the metal-rich interior expands adia- 
batically into the cavities created by the shock. 



Distribution of Metals 

The dispersal of metals by the first SN explosions 
transformed the IGM from a simple primordial gas to a 
highly complex medium in terms of chemistry and cool- 
ing, which ultimately enabled the formation of the first 
low-mass stars. However, this transition required at least 
a Hubble time, since the presence of metals became im- 
portant only after the SN remnant had stalled and the en- 
riched gas re-collapsed to high densities [42]. Further- 
more, the metal distribution was highly anisotropic, as 
the post-shock gas expanded into the voids in the shape 
of an 'hour-glass', with a maximum extent similar to the 
final mass-weighted mean shock radius (see Figure 11) 
[42]. 

To efficiently mix the metals with all components of 
the swept-up gas, a DM halo of at least Myir ~ 10^ Mq 
had to be assembled [42], and with an initial yield of 0. 1, 
the average metalUcity of such a system would accumu- 
late to Z ~ 1Q~^-^ Zq, well above any critical metallic- 
ity [15, 20, 82]. Thus, if energetic SNe were a common 
fate for the first stars, they would have deposited metals 
on large scales before massive galaxies formed and out- 
flows were suppressed. Hints to such ubiquitous metal 
enrichment have been found in the low column density 
Lya forest [3, 87, 88], and in dwarf spheroidal satellites 
of the Milky Way [47]. 



OBSERVATIONAL CONSTRAINTS 



An increasing number of direct and indirect observa- 
tions have become feasible during the last few years, 
due to rapid progress in observational methods and tech- 
niques. The most prominent among these concern the 
CMB optical depth to Thomson scattering, the near-IR 
background, high-redshift GRBs, and the possibility of 
scrutinizing the nature of the first stars by metals found 
in the oldest Galactic halo stars, dubbed 'stellar arche- 
ology'. We here briefly discuss these promising obser- 
vational avenues, and elaborate on their implications for 
the first stars and galaxies. 



Optical Depth to Thomson Scattering 

One of the most important constraints on early star for- 
mation is the CMB optical depth to Thomson scattering, 
recently revised to T ~ 0.09 ± 0.03 after the Wilkinson 
Microwave Anisotropy Probe (WMAP) three-year data 
release [89]. Combined with the absence of the Gunn- 
Peterson trough in the spectra of high-redshift quasars, 
this measurement provides an integral constraint on the 
total ionizing photon production at z > 6 [6, 102]. 

In recent work, the contribution of the first stars to the 
ionizing photon budget was determined by semianalytic 
star formation rates based on the Sheth-Tormen formal- 
ism [41]. It was found that a top-heavy primordial popu- 
lation (consisting of Pop III. 1 stars) must be terminated 
fairly rapidly in order to not overproduce the observed 
optical depth (see Figure 12). Such efficient feedback is 
possible via photoionization heating, which delays and 
possibly changes the mode of star formation in neigh- 
bouring minihalos, depending on the state of the mini- 
halo collapse [4, 63, 93, 99, 105]. 



Near-Infrared Background 

Another observation with the potential to yield valu- 
able information on the nature of the first ionizing 
sources is the detection of an excess in the near-IR back- 
ground of the order of 1 nW m^^ sr^' in the wave- 
length regime 1 — 2 ;Um [24]. If a fraction of this light 
is contributed by early stars, it corresponds to redshifted 
Lya photons emitted between z — 1 and z ~ 15. Ap- 
plying a semianalytic model of primordial star forma- 
tion and deriving the total luminosity, it has been found 
that Pop III. 1 and Pop III. 2 combined must contribute at 
a level less than 10^^ nW rcT^ sr^^ in order to avoid 
premature reionization [41]. To determine what fraction 
is contributed by metal-free stars, one must analyze the 
fluctuation power of the IR background at scales small 
enough to capture their strong clustering [29, 52, 62, 79]. 




FIGURE 12. Ionization histories (Panel (a)), and implied 
optical depths (Panel (b)) for a weak (dotted Une), intermediate 
(dashed line), and strong (dot-dashed line) Pop III.l mode. For 
comparison, Panel (b) includes the range of allowed values for 
the optical depth (long-dashed lines). 



High-redshift Gamma-ray Bursts 

GRBs are beUeved to originate in compact remnants 
(neutron stars or black holes) of massive stars, and their 
high luminosities make them detectable out to the edge 
of the visible universe [27, 56]. GRBs offer the opportu- 
nity to detect the most distant (and hence earliest) pop- 
ulation of massive stars. In the hierarchical assembly 
process of DM halos, the first galaxies should have had 
lower masses (and lower stellar luminosities) than their 
low-redshift counterparts. Consequently, the characteris- 
tic luminosity of galaxies or quasars is expected to de- 
cUne with increasing redshift. GRB afterglows, which al- 
ready produce a peak flux comparable to that of quasars 
or starburst galaxies at z ~ 1 — 2, are therefore expected 
to outshine any competing source at the highest red- 
shifts. As the electromagnetically-brightest explosions 
in the universe, GRBs should be detectable out to red- 
shifts z > 10 [27, 56]. High-redshift GRBs can be iden- 
tified through IR photometry, based on the Lya break 
induced by absorption of their spectrum at wavelengths 
below 1.216 /im [(1 +z)/10]. Follow-up spectroscopy of 
high-redshift candidates can then be performed on a IO- 
meter-class telescope. Recently, the ongoing Swift mis- 
sion has detected a GRB originating at z — 6.3 [40, 44], 
thus demonstrating the viability of GRBs as probes of the 
early universe. 



Although the nature of the central engine that pow- 
ers the relativistic jets of GRBs is still unknown, recent 
evidence indicates that long-duration GRBs trace the for- 
mation of massive stars [11, 12, 68, 95, 100], and in par- 
ticular that long-duration GRBs are associated with Type 
Ib/c SNe [91]. Since the first stars in the universe are 
predicted to be predominantly massive [1, 14, 17], their 
death might give rise to large numbers of GRBs at high 
redshifts. In contrast to quasars of comparable bright- 
ness, GRB afterglows are short-lived and release 10 or- 
ders of magnitude less energy into the surrounding IGM. 
Beyond the scale of their host galaxy, they have a negligi- 
ble effect on their cosmological environment. However, 
the feedback from a single GRB or SN on the gas con- 
fined within the first galaxies could be dramatic, since 
the binding energy of most galaxies at z > 10 is lower 
than 10^^ ergs [7]. Consequently, they are ideal probes of 
the IGM during the reionization epoch. Their rest-frame 
UV spectra can be used to probe the ionization state of 
the IGM through the spectral shape of the Gunn-Peterson 
(Lya) absorption trough, or its metal enrichment history 
through the intersection of enriched bubbles of SN ejecta 
from early galaxies [38]. 

What is the signature of GRBs that originates in metal- 
free. Pop III progenitors? Simply knowing that a given 
GRB came from a high redshift is not sufficient to reach a 
definite conclusion as to the nature of the progenitor For 
example, the currently highest-redshift GRB at z ~ 6.3 
clearly did not originate from a Pop 111 progenitor, given 
the inferred level of metal enrichment in the host sys- 
tem of a few percent solar [25]. Pregalactic metal en- 
richment was Ukely quite inhomogeneous, and we ex- 
pect normal Pop 1 and 11 stars to exist in galaxies that 
were already metal-enriched at these high redshifts [22]. 
Pop III and Pop I/II star formation is thus predicted to 
have occurred concurrently at z > 5. How is the predicted 
high mass-scale for Pop 111 stars reflected in the obser- 
vational signature of the resulting GRBs? Preliminary 
results indicate that circumburst densities are systemat- 
ically higher in Pop HI environments. GRB afterglows 
will then be much brighter than for conventional GRBs. 
In addition, due to the systematically increased progen- 
itor masses, the Pop III distribution may be biased to- 
wards long-duration events. Figure 13 leads to the robust 
expectation that ~ 10% of all Swift bursts should origi- 
nate at z > 5. This prediction is based on the contribu- 
tion from Pop 1/11 stars which are known to exist even at 
these high redshifts. Additional GRBs could be triggered 
by Pop III stars, with a highly uncertain efficiency. 



Stellar Archeology 

The discovery of extremely metal-poor stars in the 
Galactic halo has made studies of the chemical com- 
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FIGURE 13. Predicted GRB rate to be observed by Swift 
[22]. Shown is the observed number of bursts per year, 
dNQ'^^/d\n{l + z), as a function of redshift. All rates are 
calculated with a constant GRB efficiency, tjgrb — 2 x 
10^'^ bursts M^'. Dotted lines: Contribution to the observed 
GRB rate from Pop I/II and Pop III for the case of weak chemi- 
cal feedback. Dashed lines: Contribution to the GRB rate from 
Pop I/n and Pop ni for the case of strong chemical feedback. 
Filled circle: GRB rate from Pop III stars if these were respon- 
sible for reionizing the universe at z ~ 17. 



position of low-mass Pop II stars powerful probes of 
the conditions in which the first low-mass stars formed. 
While it is widely accepted that metals are required for 
the formation of low-mass stars, two general classes of 
competing models for the Pop III - Pop II transition 
are discussed in the literature: (i) atomic fine-structure 
line cooling [20, 78]; and (ii) dust-induced fragmenta- 
tion [82]. Within the fine-structure model, C II and O I 
have been suggested as main coolants [20], such that 
low-mass star formation can occur in gas that is enriched 
beyond critical abundances of [C/H]grit — —3 .5 ± 0. 1 and 
[0/H](.rit ~ — 3 ± 0.2. The dust-cooling model, on the 
other hand, predicts critical abundances that are typically 
smaller by a factor of 10 — 100. 

Based on the theory of atomic line cooling [20], a 
new function, the 'transition discriminant' has been in- 
troduced: 

i^trans = logjo (lO^C/H] + 0.3 X 10^^^"^) , (1) 

such that low-mass star formation requires Dtrans > 
Dtransjcrit — —3.5 ±0.2 [37]. Figure 14 shows values of 
Otrans for a large number of the most metal-poor stars 
available in the Uterature. While theories based on dust 
cooling can be pushed to acconmiodate the lack of stars 
with Dtrans < Otrans,crit> it appears that the atomic-cooling 
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FIGURE 14. Transition discriminant, f trans, for metal-poor 
stars collected from the literature as a function of [Fe/H]. Top 
panel: Galactic halo stars. Bottom panel: Stars in dSph galaxies 
and globular clusters. G indicates giants, SO subgiants. The 
critical limit is marked with a dashed line. The dotted lines refer 
to the uncertainty. The detailed references for the various data 
sets can be found in [37]. 



theory for the Pop III - Pop II transition naturally ex- 
plains the existing data on metal-poor stars. Future sur- 
veys of Galactic halo stars will allow to further populate 
plots such as Figure 14, and will provide valuable insight 
into the conditions of the early universe in which the first 
low-mass stars formed. 



OUTLOOK 

Understanding the formation of the first galaxies marks 
the frontier of high-redshift structure formation. It is cru- 
cial to predict their properties in order to develop the op- 
timal search and survey strategies for the JWST. Whereas 
ab-initio simulations of the very first stars can be carried 
out from first principles, and with virtually no free pa- 
rameters, one faces a much more daunting challenge with 
the first galaxies. Now, the previous history of star forma- 
tion has to be considered, leading to enhanced complex- 



ity in the assembly of the first galaxies. One by one, all 
the complex astrophysical processes that play a role in 
more recent galaxy formation appear back on the scene. 
Among them are external radiation fields, comprising 
UV and X-ray photons, and possibly cosmic rays pro- 
duced in the wake of the first SNe [90]. There will be 
metal-enriched pockets of gas which could be pervaded 
by dynamically non-negligible magnetic fields, together 
with turbulent velocity fields built up during the virializa- 
tion process. However, the goal of making useful predic- 
tions for the first galaxies is now clearly drawing within 
reach, and the pace of progress is likely to be rapid. 
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